Finite-size effects in strongly interacting Rydberg gases 
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The scaling of the number of Rydberg excitations in a laser-driven cloud of atoms with the 
interaction strength is found to be affected by the finite size of the system. The scaling predicted by 
a theoretical model is compared with results extracted from a numerical many-body simulation. We 
find that the numerically obtained scaling exponent in general does not agree with the analytical 
prediction. By individually testing the assumptions leading to the theoretical prediction using the 
results from the numerical analysis, we identify the origin of the deviations, and explain it as arising 
from the finite size of the system. Furthermore, finite-size effects in the pair correlation function 
g^'^-' are predicted. Finally, in larger ensembles, we find that the theoretical predictions and the 
numerical results agree, provided that the system is sufficiently homogeneous. 
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I. INTRODUCTION 

During the past few years, a growing amount of re- 
search, both in experiment and theory, has been aim- 
ing at a detailed understanding of ultracold gases of 
atoms excited to states of high principal quantum num- 
ber. Such Rydberg atoms have extreme properties, of 
which most important are their long-range dipolar inter- 
actions. These imply the so-called dipole-blockade pre- 
dicted in [H, [1], a suppression of further Rydberg exci- 
tation in the vicinity of one such excited atom, which 
is caused by an interaction-induced level shift. First ex- 
perimental detection of the dipole blockade was reported 
in Refs. Q and [1]. Subsequently, further effects such 
as collective oscillations and supcrradiance, to name a 
few, were found At sufficiently low temperatures, 
the atomic motion and collisions can be neglected on 
the timescale of the electronic dynamics, allowing for the 
'frozen-gas' approximation. These properties, combined 
with high control over the atomic positions provided by 
optical trapping techniques, make up an ideal model sys- 
tem for strongly correlated ensembles of atoms. Various 
applications in quantum information processing Q and 
quantum simulation Q are being pursued. With recent 
progress, also spatially resolved measurements are com- 
ing within reach @-[l3- 

As a particular consequence of the many-body correla- 
tions, the ground state of a strongly interacting and laser- 
driven Rydberg ensemble can undergo a second-order 
quantum phase transition [ll| . At the critical point, de- 
termined by resonant laser driving, the system can be 
characterized by a single parameter a which depends on 
the laser Rabi frequency fl, the ensemble density and 
the interaction strength. The number of excited Ryd- 
berg atoms A^Ryd was predicted to scale as A^Ryd ^ a", 
with an exponent in the critical region a <C 1. This re- 
sult was found to agree well with numerical results from 
mean-field calculations and many-body simulations 
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In this paper, we revisit this scaling of the amount 



of Rydberg excitation with the interaction strength and 
show that it can be modified by finite-size effects in 
smaller ensembles of Rydberg atoms. For this, we com- 
pare the scaling prediction obtained from a theoretical 
model with results extracted from a numerical many- 
body simulation. We find algebraic scaling in both the 
model and the simulation. But the numerically obtained 
scaling exponent v in general does not agree with the an- 
alytical prediction. By testing the assumptions entering 
the argument leading to the prediction we identify the 
deviations to arise from the finite size of the system. To 
be precise, the estimation of the blockade radius based 
on the comparison of the two energy scales present in the 
system turns out not to be biased by the finite size of the 
system. On the other hand, relating the size of the block- 
ade spheres to the steady-state excited fraction is found 
to be non-trivial in finite-size systems. In this context we 
discuss the validity of the so called 'super-atom' picture 
in finite systems. We also find finite-size effects in the 
pair correlation function g^'^\ We finally analyze larger 
ensembles, imposing periodic boundary conditions. We 
find that, provided that the ensemble is sufficiently homo- 
geneous, finite-size effects disappear and the numerically 
obtained scaling exponent agrees well with the predicted 
one. We point out that, while finite-size effects appear 
in the scaling exponent of the bulk number of Rydberg 
atoms, their interpretation crucially relies on a spatially 
resolved analysis of the excitation dynamics. 

Our study is not least motivated by the fact that a 
number of recent experiments approach length scales for 
which finite size effects are expected to become rele- 
vant [T3 - [l6j . In particular, the system sizes considered 
in the following are well within reach of experimental re- 
alizations of quasi ID Rydberg ensembles |15| . 

After introducing the model and numerical methods in 
Sect, [n] we analyze, in Sect. IIIIl the dynamical evolution 
of the ensemble under the influence of the laser radia- 
tion. In Sect. IIV Al we study bulk properties of the en- 
semble and show that the numerical simulations predict 
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scaling parameters different from tfie analytical model. 
To identify the origin of this deviation, we compute spa- 
tially resolved observables in Sect. lIVBl In Sect. lIVCl we 
compare to systems with periodic boundary conditions. 
Section llV Dl deals with finite-size effects in the pair cor- 
relation function. We summarize our results in Sect. |Vl 



II. MODEL AND METHODS 
A. Hamiltonian 

Our model system is a cloud of two-level atoms con- 
sisting of a ground state \g) and a Rydberg state |r). The 
atoms reside at fixed positions throughout the simulation 
time. In most experiments, the atoms are excited from 
\g) to |r) via a two- photon transition coupling to an in- 
termediate level |m). However, if the intermediate level 
is far detuned it can be adiabatically eliminated [l7j . In 
the present work we assume that the intensity and wave- 
length of the laser is temporally and spatially constant. 

The Hamiltonian describing such a cloud of N in- 
teracting, laser-driven two-level atoms can be written 

as U,um 



B. Time evolution 
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Here, Hi^ is the one-body Hamiltonian of N non- 
interacting atoms, written in rotating wave approxima- 
tion and in units where ft = 1 as 
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where we used the notation s^*^ = \a)^ for the internal 
projector of the ith atom. H^. depends on the detuning 
A = Ai + A2 = ujii + LO12 — LOgrn — i-Ornr between the laser 
frequencies luh and the energy differences oJq,^ = — E^ 
between the ground and intermediate, and the interme- 
diate and Rydberg states. H^i contains the two-photon 
Rabi-frequency f2 that depends on the laser intensities 
and the dipole matrix elements of the transitions. We 
have assumed to be real. 

The long-range interactions between two Rydberg 
atoms are accounted for by the van der Waals term 
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In our simulations, we set fl — 1, i.e., we measure 
energy in units of ftfi and time in units of 1/0. We keep, 
however, Cg as a dimensionful parameter and give lengths 
in micrometers, keeping in mind that increasing Cg is 
equivalent to decreasing the volume at a fixed particle 
number, i.e., increasing the density. 



We numerically solve the Schrodinger equation 
z||vl/)=i/|*) 



(4) 



for a set of A'' two-level atoms at positions and given 
parameters Cg and A, starting from an initial config- 
uration where all atoms are in the ground state. A 
variety of approaches has been suggested to deal with 
the large dimension of Hilbert space of the above sys- 
tem which scales exponentially as d ^ 2^ . These in- 
clude techniques relying on an expansion in the order 
of correlations [3, [ll|, [l^, [l^l, rate-equation methods 
[21I - I23 I , and man y-b ody simulations on truncated state 
spaces pJl [iSl |25| - |35| . In the following analysis, we trun- 
cate the Hilbert space, making use of the fact that states 
1$) with high diagonal elements (<I>|i/int + \^) are 
less likely to become populated. We note that for infi- 
nite interaction strength, starting from the ground state 
|0) = \gg . . .g), only singly excited states are accessible. 
Higher excited states can be discarded because they have 
infinite interaction energy. The resulting dynamics is a 
Rabi oscillation between the ground state |0) and the 
symmetrized singly excited state |s) = -/V"^/^^-^-^ 

with collective Rabi frequency fic — VND. [2^. Col- 
lective Rabi oscillations of two atoms could already be 
observed in experiment [s^ . 

Taking into account the dipole-dipole interactions 
which inhomogeneously shift the transition frequency, 
two atoms cannot be excited simultaneously within the 
blockade volume Vh- All atoms in this volume share 
one excitation and their dynamics is governed by the 
collective Rabi-frequency ^/NlVl. The blockade volume 
is found by comparing the relevant energies CqR^^ ~ 
\/Nhfl where the blockade radius Ri, is the minimal dis- 
tance that two simultaneously excited atoms can have. In 
a d-dimensional gas of density p (atoms per unit length 
to the power d) we have Nj, = pVb oc pR'^- These consid- 
erations lead to [1, 0, [ll| 



Rb oc 



2/(d+12) 



(5) 



Note that in deriving R^, we neglected the contribution 
of Ha to the diagonal element. The above estimate is 
therefore only valid near resonance, i.e., for A ^ il. 

We use Eq. ^ to truncate the Hilbert space by ne- 
glecting all states in which two excited atoms have a dis- 
tance smaller than the cutoff distance < Rb- This 
is achieved by choosing the proportionality constant in 
Eq. ([5|) appropriately. Other ways to reduce the number 
of basis states are to put an upper bound 

^max the 

number of excitations or to exclude states with an in- 
teraction energy higher than some threshold S;™^''. We 
have checked that our numerical results are converged 
with respect to all truncation criteria, i.e., to the chosen 
values of the cutoff parameters. 
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To implement the truncation we construct, from the 
canonical set of product states |Q;ia2 • ■ • Q^w) = |cti) (8) 
|a2) CS) . . . lajy) with € {g, r}, basis states of the form 
|m; ii, i2 . . . im), i-e., states with m atoms ii, 12 ■ ■ ■ im hi 
the Rydberg state and all other atoms in the ground 
state. The detuning and interaction parts of the Hamil- 
tonian are diagonal in this basis, while the Rabi terms 
connect states that differ in the internal level of one atom 
only. Therefore the Hamiltonian matrix can be divided 
into sub blocks consisting of states containing the same 
number m of excited atoms. The Hamiltonian has a tridi- 
agonal structure in these blocks since only states with m 
differing by one are connected by the laser Hamiltonian. 
The blocks on the diagonal are diagonal in the above 
basis. This structure allows for a convenient truncation 
and construction of the resulting evolution matrices. Af- 
ter the truncation, we denote the n remaining basis states 
of form |m; ii, 12 . . . im) as |$fc) with k G {1, . . . , n}. 

While the performance of our code depends on the 
number of atoms it is crucially affected by the number 
of excitations that fit into the trap volume. The major 
numerical limitation is computation time. If we want to 
limit the time needed for the propagation of a single re- 
alization into the steady state to less than one hour on 
a state-of-the-art desktop processor (2.33 GHz), we are 
limited to state spaces of about four million states. Given 
this constraint, if we consider an atom number of about 
100, the trap size needs to be limited such that it can 
contain at most six-fold excitations. For 60 atoms we are 
limited to about nine excitations. For such system sizes 
the coherent time evolution of the system does not fully 
converge to a steady state, with small oscillations persist- 
ing in the long-time limit. Moreover, for the calculation 
of spatially dependent observables such as the distribu- 
tion of the Rydberg excitations, a single run with about 
100 atoms does not provide sufficient statistics. These 
problems can be overcome by artificially increasing the 
particle number by Monte Carlo sampling. This simply 
means that we calculate the time evolution of the sys- 
tem for several realizations (typically 1000 to 5000) with 
randomly chosen atom positions and average over the 
resulting observables. This method was used before in 
Refs. [13, 113, and can be justified for Bose-Einstein 
condensates as well as for thermal clouds as explained in 
Ref. [H. 

Our model exclusively contains the unitary time evo- 
lution, i.e., spontaneous emission, dephasing or laser 
linewidths are not explicitly accounted for. A proper 
treatment of such incoherent processes would require 
methods such as a master equation or a quantum 
Monte Carlo simulation of the incoherent jumps, which 
would considerably increase the computational complex- 
ity. From a physical point of view, this approxima- 
tion can be justified if the experimental time scales are 
short compared to the life time of highly excited Rydberg 
states, and if stabilized laser systems with low linewidth 
are used. 



C. Observables 

For a single realization of our system, all information 
is contained in the coefficient vector c with elements Ck, 
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with basis states as defined in Sec. IHBI For each 
realization the observables are calculated from the state 
vector at a given time of interest and then averaged by 
Monte Carlo sampling. The mean number of Rydberg 
atoms is 



n 



k=l 



(7) 



where to'^*^) is the number of excited atoms in state \^k) 
The excitation probability of the ith atom is given by 



k=l 



ICfcl 



(8) 



where a,-'^^ is the internal state of atom i in state \^k), 
and the Kronecker delta selects the coefficients of the 
states in which the zth atom is excited. 

The spatial distribution of the Rydberg excitations 
A^ryd(r) is obtained by dividing the trap volume into 
small cells and summing up the excitation probabilities of 
all atoms that lie in such a cell, normalized by their num- 
ber. Most important in this work are spatial correlations 
quantified by the pair correlation function ^'•^•'(r, r') [40| . 
This function is a measure for the conditioned probabil- 
ity of having an excitation at r if there is already one 
at r'. Since the interaction potential in our model only 
depends on the mutual distance of two atoms, it is valid 
to assume that also g^"^^ only depends on r = |r — r'|. 
This has been verified numerically. The pair correlation 
of two particles is defined as (4^ 
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Discretizing space, we define 



E(r, Ar) 



(2) 



1 



(10) 



where 'Y^jj''/^^^ denotes the sum over all pairs with mutual 
distance |ri — rj| lying within the interval [r, r-f-Ar]. This 
means that we sum up the correlations of all these pairs 
divided by their number. Note that some authors use 
the alternative definition g'-'^' 



to our definition g^'^^ (r) 
atoms. 
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1 corresponds to uncorrclatcd 
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FIG. 1. Time evolution of the number of Rydberg excitations 
in an ensemble of = 50 atoms enclosed in a cylindrical trap 
with radius r = 1 fj,m and length L = 10 /im. Further param- 
eters are: A — 0, Ce ~ 900 Q fim^. (a) Single realization. The 
remaining oscillations at long times reflect the fluctuations of 
the positions of the Rydberg atoms, (b) Average over 1000 
Monte Carlo samples. 



FIG. 2. (Color online) Number A'Ryd of Rydberg excita- 
tions in a one dimensional gas. (a) A''Ryd as a function of 
the line density. Parameters: L = 10 /im, Ce ~ QOOfi/im®, 
A = 0. (b) Number of excitations as a function of the 
van der Waals interaction constant, on a double logarithmic 
scale. Crosses: numerical simulation. Solid line: power-law 
flt A^Ryd ~ Cg""-^^**. Parameters: L = 15 /im, iV = 45, A = 
of A'' = 45 atoms. 



III. DYNAMICAL PROPERTIES 

A single atom or a gas of non-interacting atoms under- 
goes simple Rabi oscillations on a time scale short com- 
pared to the dephasing time. If the laser is detuned, only 
a part of the population oscillates between the ground 
and excited states, while the rest of the population re- 
mains in the ground state. The oscillation frequency is 
modified to -y/A^ + Jl^. In the case of infinitely strong 
interactions the system reduces to a two-level system 
showing collective Rabi oscillations between the ground 
state and the collective singly excited state as discussed 
in Sec. ini 

In the intermediate case, in which several atoms can 
be excited simultaneously, but interactions are non-zero, 
the dynamics is more complicated. In the example shown 
in Fig. [1] we use a cylindrical trap of length 10 fj,m and 
diameter 2 fim. We observe that the Rydberg population 
of the cloud, initially being in the ground state, shows a 
saturation behavior. However, even for long times some 
oscillations remain, cf. Fig. [IJa). The strength of these 
oscillations depends on the number of excitations and on 
the number of atoms in the trap. This saturation can be 
interpreted in two ways: On the one hand, the mean-field 
picture predicts that, due to the disorder induced by the 
randomness of the atom positions, interaction shifts S dif- 
fer from atom to atom, leading to different Rabi frequen- 
cies y/S^ + D^ for the individual atoms. Consequently, 
the oscillations dephase over time and lead to a satura- 
tion of the overall excitation. Alternatively, one uses the 
so-called "super-atom" picture: Any excited atom in the 
ensemble blocks the excitation of the surrounding atoms 
which leads to collective excitation of the atoms within 
this blockade sphere. Again due to the random positions 
of the atoms, the number of atoms per blockade sphere 
and thus the collective Rabi frequency varies across the 
ensemble. The occurrence of different collective Rabi fre- 
quencies leads to a dephasing and thus to saturation. 



Averaging over many realizations shows that the mean 
number of Rydberg excitations saturates fully after some 
time, cf. Fig. [Ub). The timescale of this saturation de- 
pends on how many Rydberg excitations are present in 
the gas, on the system size and also the detuning. Note 
that for very small traps that confine the gas to near 
one blockade volume, ongoing strong oscillations are ob- 
served even after Monte Carlo averaging. The reason is 
that in this case, only one or a few different collective 
Rabi frequencies are possible, depending on the number 
of atoms in the blockade spheres. In the example shown 
in Fig.[IJb) we find that the excitation number converges 
to about A^ryd = 2.5 Rydberg excitations. This means 
that Nb = N/N^yd = 20 atoms share one excitation 
and we expect to observe a collective Rabi-frequency of 
flc/i^ = y/Nb = 4.47. Measuring the period of the initial 
oscillations in Fig. [TJb) wc obtain ilc = 2tt/T = 4.49 f2, 
consistent with the super-atom picture. 



IV. FINITE-SIZE EFFECTS 

A. Signatures in bulk properties 

In this section we study the properties of the system 
after is has saturated to its steady state. For sufficiently 
high atomic density the number of Rydberg excitations 
does not scale linearly with the number of atoms, as it 
would be expected for a dilute non-interacting gas. In- 
stead, a reduced Rydberg population is observed as is 
shown in Fig. [2[a). This is caused by the well-known 
dipole blockade l3| inhibiting the simultaneous excitation 
of two nearby atoms by means of a large interaction shift. 
This effect is the more pronounced the stronger the in- 
teractions arc. Thus A'ryd decreases as a function of the 
interaction parameter Cg, as is illustrated in Fig. [5l^b). 
Assuming that the blockade radius is given by Eq. ([5]), 



FIG. 3. (Color online) Pair correlation function defined in 
Eq. (|f op for three different interaction strengths Ce after t = 
200 SI evolution time. One-dimensional trap of length L — 
15 /im holding A'^ = 45 atoms exposed to a resonant laser. 



FIG. 4. (Color online) Pair correlation function (jlOp as a 
function of (a) the atom line density, with Ce = QOOflfim^, 
L — 10 fim; (b) the interaction constant Ce, with line density 
p = 3//im (L — 15 nm, N = 45). Dashed lines indicate the 
estimates for Rb as in Eq. (|12p. 



one finds that (ill 
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From this, for a one-dimensional gas (d = 1), one ob- 
tains an algebraic scaling of iVi-yd with Cq , with exponent 
1/ = -2/13 « -0.1538. Fig. |I[b) confirms that the ex- 
cited fraction scales as a negative power of Cg, however, 
with a fitted exponent i/fn = — 0.1283 ± 0.0006 (errors we 
give on fitted values just refer to the confidence interval 
of the fit). This is obtained for a one-dimensional trap of 
length 15 /im containing 45 atoms. At first sight this rc- 
sult appears to be in contrast to the calculations reported 
in [11| where good agreement of numerical simulations 
with the analytically predicted value was reported. We 
will analyze this discrepancy in more detail in the fol- 
lowing section and show that the deviation from the the- 
oretical prediction is due to finite-size effects arising in 
our relatively small trap volume. For now we summarize 
that the excited fraction follows an algebraic scaling over 
a wide range of interaction strengths also in finite sys- 
tems, while the scaling exponent deviates from the one 
predicted for an unbounded homogeneous gas. 



B. Origin of the finite-size effects 

In the previous subsection we reported the numerically 
obtained power-law scaling of the number of Rydberg 
excitations A'^ryd ~ Cq. The value for v did not agree, 
however, with the value predicted and confirmed previ- 
ously for a homogeneous gas. In the following, we trace 
back the origin of this discrepancy to finite-size effects. 
For this, we individually examine the two assumptions 
which lead to the prediction of the scaling parameter 
in Eq. ([TT|). The first one is the estimate of the block- 
ade radius, Eq. ([5]). The second one is the assumptions 
that the number of Rydberg excited atoms can be ob- 
tained from an estimate based on the super-atom picture, 
A'^ryd oc mmax = N/Nb = L/Rb- Notc that the effects 
studied here are physical finite-size effects in the sense 



that they originate from the finite trap volume rather 
than from computational limitations. Therefore, they 
could possibly be observed in experiments. 



1. Estimate for the blockade radius 

In order to analyze the estimate for the blockade 
radius given in Eq. ([5|) we consider pair correlations. 
Fig. [3] shows the pair correlation function pUj) of a one- 
dimensional gas in its thermalized state. Pairs with small 
mutual distance arc never excited simultaneously, hence 
the g'-^'-function is zero for small r. Note that the numer- 
ical blockade radius rf, was chosen well below the onset 
of non-zero correlations to ensure that the state space 
truncation does not affect the outcome of the simula- 
tions. The sharp peak which emerges just outside the 
blockaded region is located precisely at 



Rb — 



2/13 



(12) 



which was our estimate of the blockade radius (Eq. ([5])) 
leading to Eq. (Hi}. In order to probe this conjecture 
more thoroughly we analyze g^^^ for varying line density 
and interaction strength. Results are shown in Fig. 2] 
The predicted blockade radius perfectly matches the first 
maximum of g^^-* in the whole parameter region. This in- 
dicates that the estimate for Rb in Eq. ([5]) are reasonable 
even in finite geometries. We remark that in these figures 
one can also observe that correlations become stronger 
as the density or the interaction strength increases. We 
could also quantitatively verify Eq. ([T^ . For this, we ex- 
tracted the blockade radius from the simulation results 
shown in Fig.|3]by defining it as the pair distance at which 
the f/'^-* -function first exceeds one. This observable shows 
a power law dependence over the whole range of param- 
eters and fitting the exponent yields perfect agreement 
with Eq. ([Ta . 
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FIG. 5. (Color online) Same setup as in Fig. [3] (a) A'ryd 
as a function of L/Rb = L{C,i/Q.^)~'^^'^'^^ . Red line: Nu- 
merical data, obtained by varying Cs- Dotted black line: 
0.5 L/i?6. Dashed black line: 0.5(L/i?i, + 1). Solid black line: 
0.5(I/cff/-R6 + !)• None of these estimates fits the numerically 
obtained dependence, (b) Spatial distribution of the Rydberg 
excitations on the string in excitations per ^ra, divided by the 
total number of excitations. 



2. Estimate for the number of Rydberg excited atoms 

Next, we consider the assumption that A^ryd ~ fnmax = 
N/Nh ~ L/Rfj. In Fig. EJa), results are shown for the 
numerically obtained A^ryd a-s a function of Cg , expressed 
in units of L/R), via Eq. The dotted line shows 

^*L/Ri,, where 7* is the average fraction to which the 
super- atoms are excited and was assumed to be 0.5. It 
can be seen that this estimate does not coincide with the 
numerically found iVj-yd- 

This failure in finite systems is due to two reasons. 
First, from Fig. [SJb) it can be seen that the distribu- 
tion of excitations on the string shows strong maxima at 
both ends. These can be understood in the mean-field 
picture. The interaction shift of atoms close to the bor- 
der is smaller than that of atoms in the center since they 
only have potentially excited neighbors to one side. Con- 
sistent with the structure of the pair correlation function, 
the side maxima have to be followed by minima due to 
the blockade effect and, towards the center of the string, 
by further maxima corresponding to the first maximum 
of g^'^\ The maxima at the edges of the string lead to 
a higher excited fraction compared with the case of an 
infinite string. This effect becomes more dominant as Cg 
increases since the number of excitations on the string 
decreases leading to an effectively smaller system. 

Secondly, in a system of length L, the assumption 
of densely packed excitations at distances Rb from each 
other docs not necessarily lead to TVryd = 'j*L/Rb. The 
number of super-atoms fitting on such a string is rather 
L/Rb + 1 due to the excitations sitting at both ends of 
the string. This would result in iVj-yd = 7*(L/i?;, + 1) 
which fits the numerical data much better but overesti- 
mates it a bit (dashed black line in Fig.[SJa)). A further 
obvious effect is connected to the finite density of the 
gas. The finite atomic line density of p = 3/xm~^ im- 
plies that it is very unlikely that two atoms are sitting 
exactly at the trap ends, so the outermost super-atoms 



are shifted towards the trap center on average by 0.5 
leading to an efficient trap length of Lcff = L — . We 
therefore show N^y^ = -j* {Lcff/ Rb + 1) (solid black line 
in Fig. EJa)) where again 7* ~ 0.5 was assumed. This 
line now reproduces the slope of the numerical data well 
but still overestimates it. In order to understand this, 
we should bear in mind that the notion of super-atoms 
sitting at the edges of the trap is not consistent with the 
ansatz we used to estimate the blockade radius. Rb is 
estimated by equating the collective Rabi frequency of 
Nb atoms to the interaction strength at a distance Rb, 
i. e., ^/NbQ = Ce/Rb- A super-atom sitting at the edge 
of the trap would now contain less atoms, leading to a 
larger blockade radius. We conclude that the super-atom 
picture is not adequate to explain the effects of the trap 
boundaries. 

In general, the validity of universal scaling laws de- 
pends on how well length scales are separated from each 
other. The blockade radius should be much smaller than 
the trap size but much larger than the intermediate par- 
ticle distance [llj. As the simulation in Fig. [S] was done 
at constant atom density and trap size, increasing L/Rb 
means that the number of atoms per blockade sphere de- 
creases. Thus by increasing L/Rb we effectively move 
towards the regime of low density, where the interatomic 
distance is of similar size as the blockade radius. In this 
regime the predicted scaling laws are not expected to 
hold. 

Fig. [5ja) shows yet another interesting property of 
small clouds. If L/Rb is smaller than one, the cloud is 
perfectly blocked and Rabi oscillates between the ground 
state and the fully symmetric singly excited state as men- 
tioned in Sect. IIIIl Thus the time average of iVryd is 0.5 
as confirmed by the simulation data. As L/Rb increases, 
according to the super-atom picture, the number of ex- 
citations should immediately jump to A^ryd = 1 as soon 
as two super-atoms fit into the trap. However, A^ryd in- 
creases smoothly, indicating that super-atoms cannot be 
viewed as hard spheres but are rather soft objects. It 
should also be mentioned that, as the trap size is in- 
creased holding density and all other parameters con- 
stant, for L/Rb > 5.5, A'ryd increases linearly with a slope 
of 0.5/ Rb confirming the picture of soft super-atoms. 

Our above findings explain why the exponent v is un- 
derestimated if finite-size effects are neglected, and why 
it is quite difficult in finite systems to find an analyt- 
ical relation between A'ryd and mmax based on simple 
assumptions. For larger systems, these finite-size effects 
become negligible relative to the bulk properties of the 
large ensemble. We thus conclude that the assumptions 
A'ryd ^ TOmax = N/Nb = L/Rb fails for small ensembles, 
leading to the discrepancy between the numerically ob- 
tained scaling parameter v and the theoretical prediction 
in Eq. (tTTI) . 
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FIG. 6. (Color online) Same setup as in Fig. [S] but with 
periodic boundary conditions, (a) Pair correlation function as 
a function of Ce- Effects of artificial self coupling due to finite 
system size are visible particularly at high Ce- (b) Crosses: 
Simulated number of Rydberg excitations. Solid line: Fit of 
algebraic decay. Only the first nine points are used for the fit. 
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FIG. 7. (Color online) Finite-size effects in the g'-^'-function. 
If the correlation length is on the order of the trap size, cor- 
relations are enhanced by the finite-size effects. Parameters 
are Ce = 900 MHz fim'^ and p = 3.33 fim'^ . 



D. Finite-size effects in the correlation function g' 



C. Large ensembles via periodic boundary 
conditions 



If periodic boundary conditions are used, the excita- 
tion density is perfectly flat. However, if strong correla- 
tions are present, large system sizes are required in or- 
der to eliminate artificial self-coupling effects. The lat- 
ter are also known as aliasing effects that are caused by 
atoms correlated with a neighboring atom twice, inside 
the string and across the string boundary. Therefore the 
string length has to be at least twice as large as the range 
of the correlations. This is illustrated in Fig. El^a). At 
large values of the interaction strength, for which the 
blockade radius is large and correlations are long range, 
the g(^''-function deviates strongly from the case with- 
out periodic boundaries. Only at small values of Cg the 
system seems to become free of artificial self-coupling ef- 
fects. Note that in contrast to the finite-size effects with- 
out periodic boundary conditions discussed above, these 
self-couplings cannot be observed in a linear ID ensem- 
ble. They could, however, be observed if instead the gas 
was arranged in a suitable ID ring trap. To avoid the 
self-coupling effects, we fitted the algebraic scaling law 
to the nine points with lowest Cg shown in Fig. HJb) 
and obtained ly « —0.146 ± 0.003. The line density is 
p = 3 fim~^ as in the fixed-boundary case. For larger 
densities, i.e., larger Nf,, the scaling exponent obtained 
from our numerical simulations approaches the predicted 
value of = —2/13. This can be understood since at 
high densities, the nearest-neighbor distances between 
ground-state atoms is much smaller than Rb while the 
trap length is much larger than Ri,. In this case, the sys- 
tem appears homogeneous, as neither the finite trap size 
nor the coarseness of the atom distribution are relevant. 



We finally analyze finite-size effects in the pair corre- 
lation function. Fig. [Tl^a) shows g^'^^r) for three differ- 
ent trap lengths and a density oi p = 3.33 /im~^. We 
notice that the deviation from g^'^\r) = 1, indicating 
spatial correlations, decrease slightly with increasing L. 
We now focus on the first maximum of the correlation 
function, c. f. Fig. [Zl^a). This maximum is the feature 
that depends on the trap size most strongly. Looking 
at this peak more closely, we find a double-peak struc- 
ture which vanishes in the limit of large very trap size. 
Splitting up the (/'■^■'-function into contributions from the 
different subspaces with definite excitation number m we 
find a very pronounced peak in the contribution of the 
two-fold excited states, which is located at the position of 
the left peak sub-peak of the first maximum in Fig. [Tj^a) . 
This structure is present for any trap length, but as the 
relative population of states with a very low number of 
excitations decreases with increasing trap size, the ef- 
fect on the structure of the total g'^^ becomes weaker 
for large traps. This explains why the height of the left 
sub-peak decreases with increasing trap size, whereas the 
right sub-peak remains constant in size. The structures 
seen on the m = 2 subspace can partly be interpreted in 
the two-atom picture as arising from different excitation 
channels, as done in Ref. [l^l for the case of A ^ SI. 
However, for the = case one expects collective effects 
to be dominant, and thus it is clear that the structure 
cannot be fully understood from the perspective of the 
two-atom problem. As expected, correlation peaks are 
generally enhanced if the system size is smaller than the 
range of the correlations. Fig. [71[b) confirms that .g^^-'(r) 
converges to a smooth function in the limit of large L. In 
the simulations with largest system sizes in Fig.[71[b), the 
gr^^-'-function saturates to 1 over a distance smaller than 
the trap length. Nevertheless, some features of the cor- 
relation function (visible in Fig. as discussed above) 
still depend on the trap length. 
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V. CONCLUSIONS 



In summary, we have studied the time evolution and 
the steady-state statistics of the number 7V,.yd of Rydberg 
excited atoms in cyhndrical and onc-dimcnsional trap 
geometries, with and without periodic boundary condi- 
tions. Wc predict finite-size effects in the scahng of A^ryd 
with the van dcr Waals coefficient Cg which quantifies 
the dipolc-dipole interaction strength between the ex- 
cited Rydberg atoms. Our numerical analysis shows that 
TVryd and Cq are still connected by an algebraic scaling, 
but with a scaling exponent v different from the analyt- 
ically predicted value well-known in the literature. The 
analytical prediction relies on two assumptions, an esti- 
mate of the blockade radius, and an estimate of the num- 
ber of excited Rydberg atoms based on the super-atom 
model. We could show that the estimate of the block- 
ade radius provides a reliable prediction of the position 
of the first maximum of the pair correlation function ob- 
tained from the numerical simulations over a wide range 
of atom densities as well as interaction strengths. In con- 
trast, the analytical estimate of the number of excited 
Rydberg atoms did not agree with our numerical results. 
This was found to be due to the fact that atoms at the 
ensemble borders can only interact with other atoms to 
one side, leading to an overall enhanced Rydberg exci- 
tation. Additionally, for a continuously varying ensem- 
ble size, the number of hard-sphere super-atoms fitting 



into the volume is not definite, but rather lies between 
L/Rb and + 1. We also predicted finite-size effects 

in the pair correlation function g^'^\ for which wc have 
shown that some of its features depend on the ensem- 
ble size even for lengths L exceeding the distance range 
over which g^'^\r) is different from one. In larger ensem- 
bles, which we have simulated using periodic boundary 
conditions, we found that the scaling exponent v agrees 
with the analytical prediction, provided that the density 
is sufficiently high. Only then, the atomic distribution 
is sufficiently homogeneous. We conclude that two con- 
ditions are required to eliminate finite-size effects: First, 
the condition ^ needs to be fulfilled, such that 
the microscopic structure of the atom distribution be- 
comes irrelevant. Second, L ^ Rf, is needed in order 
to eliminate the effects of the trap boundaries. Wc have 
shown that for trap sizes oi L/Rb < 6 finite size effects 
significantly modify the dependence of the number of Ry- 
dberg excitations on the interaction strength. For typi- 
cal experimental parameters {Rt = 5 /zm to 10 /um) this 
corresponds to L = 30 /xm to 60um which is an experi- 
mentally relevant length scale |12l-[l6|. 
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